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Abstract 

We use high-redshift type la supernova and compact radio source data in order to test the infrared 
(IR) fixed point model of the late Universe which was proposed recently. It describes a cosmology 
with a time dependent cosmological constant and Newton constant whose dynamics arises from an 
underlying renormalization group flow near an IR-attractive fixed point. Without any finetuning 
or quintessence field it yields flu = = 1/2. Its characteristic i 4//3 - dependence of the scale 
factor leads to a distance-redshift relation whose predictions are compared both to the supernova 
and to the radio source data. According to the x 2 test, the fixed point model reproduces the data 
at least as well as the best-fit (Friedmann-Robertson- Walker) standard cosmology. Furthermore, 
we extend the original fixed point model by assuming that the fixed point epoch is preceded by an 
era with constant G and A. By means of a Monte Carlo simulation we show that the data expected 
from the forthcoming SNAP satellite mission could detect the transition to the fixed point regime 
provided it took place at a redshift of less than about 0.5. 
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I. INTRODUCTION 



.„ recent yea,, renovation group (RG) techniques have been e*ten 8 ive. y applied 
in cosmology, for instance in a classical averaging scenario [1| or in discussing the late- 
time behavior of the classical Einstein equations |2j. As for possible quantum gravitational 
effects, an exact functional RG equation |3| for Quantum Einstein Gravity (QEG) has been 
introduced 0, and, within certain approximations (truncations of "theory space"), its 
predictions for the scale dependence of Newton's constant G and the cosmological constant 
A have been worked out y, IMJp . 

In particular, it was found y, 10, Q] that the ultraviolet (UV) behavior of QEG is 
likely to be governed by a non-Gaussian RG fixed point which would render the theory 
nonperturbatively renormalizable. If so, QEG is mathematically consistent and predictive 
at arbitrarily short distances and, in cosmology, at times arbitrarily close to the initial 
singularity. In ref. 3] we investigated how the standard Friedmann-Robertson- Walker 
(FRW) cosmology gets modified when the scale dependence of G and A is taken into account. 
We "RG improved" Einstein's equation by replacing G — > G(k), A — > A(k), where k is the 
running mass scale which may be identified with the inverse of the cosmological time in a 
homogeneous and isotropic Universe, k oc 1/t. In this manner the RG running gives rise to a 
dynamically evolving, time dependent G and A Q]. The improvement of Einstein's equation 
can be based upon any RG trajectory k (—> (G(k), A(k)) obtained as an (approximate) 
solution to the exact RG equation of QEG. In particular, if the non-Gaussian UV fixed 
point predicted by all known solutions does indeed exist, G{k) and A(k) scale in a very 
simple manner as k — > oo. In fact, the dimensionless Newton constant g(k) = k 2 G(k) and 
cosmological constant X(k) = A{k)/k 2 are attracted towards their fixed point values g* and 
A*, respectively, and therefore, in the fixed point regime, 

G(k) = | (la) 
A(k) = Kk 2 (lb) 

In ref. [3] the improved cosmology resulting from the trajectory (JJJ was analyzed in detail. 
As the fixed point (<?*, A*) is approached for k — > oo, it applies to the very early Universe 
(t — > 0). It turned out that the cosmology of this "Planck era" is described by an essentially 
unique attractor solution to the cosmological evolution equations which might provide a 
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solution to the horizon and flatness problems of standard cosmology without invoking an 
inflationary era 1 . 

One of the remarkable properties of the attractor solution is that it dynamically adjusts 
the vacuum energy density p\ = A/87rG so as to equal precisely the matter energy density. 
In units of the critical density, the spatially flat solution has 

n A (t) = n M (t) = x - (2) 

at any time during the fixed point era. This adjustment mechanism has led to the speculation 
[lfl that not only the very early, but also the very late history of the Universe is described 
by a RG-improved Einstein equation based upon the trajectory (0). In this scenario one 
postulates that, in addition to the UV-attractive fixed point discussed so far, there exists a 
second fixed point (g* R , A* R ) in (g, A)-space towards which every trajectory k i— > (g(k), X(k)) 
within a certain basin of attraction is attracted for k — > 0. On very large scales the cutoff 
identification k oc 1/t should still be correct so that for t — > oo the IR fixed point determines 
the asymptotically late cosmology. The IR fixed point hypothesis implies that for k —>■ 
or i ^ oo, too, G and A evolve accordingly to (In the sequel we shall omit the 

superscripts "IR" from g* and A* but it should be kept in mind that the values of g* 
and A* are different at the two fixed points.) The main motivation for this hypothesis 
comes from the fact that, provided the Universe is spatially flat and has entered the fixed 
point regime already, the densities Q\ and f^M are unambiguously predicted to assume the 
value 0.5. This is intriguingly close to the values Q\ ~ 0.7, ~ 0.3 favored by recent 
observations when interpreted within standard FRW cosmology. Hence the IR fixed point 
might provide a natural solution to the "cosmic coincidence problem" jl^l which does not 
require quintessence field. 

As for their actual verification within QEG, the UV and the IR fixed points have a rather 
different status for the time being. While the existence of the UV fixed point has actually 
been demonstrated within certain approximations whose reliability was tested in detail by 
means of rather involved calculations ly, Ml , to date there is no comparable evidence for the 

n ri 

IR fixed point. (See refs. |8j, 11 81] . however.) The technical reason is that it is extremely 
difficult to follow the RG flow from the UV all the way down to the IR. Therefore a k — ► 



1 A similar RG improvement of black hole spacetimes can be found in refs. 0,0]- 
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RG-running of the form ((TJ has the status of a conjecture for the time being; it is motivated 
by its remarkable phenomenological success. We also emphasize that the RG-running of G 
and A should not necessarily be thought of as arising from quantum effects; in principle it 
could also result from a purely classical averaging (coarse graining) Q . 

T he g eneral properties of the IR fixed point (IRFP) cosmology were first discussed in 



ref. 



16 1 and will be reviewed briefly in subsection III Al In ref. [l9j] the theory of small 



perturbations about the IRFP cosmology has been developed and its consequences for the 
cosmological structure formation were discussed. 

In the present paper we are going to test the IRFP cosmology by confronting its pre- 
dictions with the currently available data on distant type la supernovae and compact (~ 
milliarcsecond) radio sources. The relevant quantities which can be compared to the data 
are the luminosity distance di(z) and the angular diameter distance (^(V), as functions of 
the redshift. In standard FRW cosmology these functions depend on Q M and fl^, and these 
parameters are fitted so as to reproduce the observations as well as possible. In the IRFP 
cosmology instead, the functions di{z) and cIa(z) contain no free parameters at all. Hence 
an acceptable fit to the data would be a rather non-trivial success of the IRFP scenario. 

In particular we apply the y 2 test to the combined dataset of 92 type la supernovae 
discovered by the Supernova Cosmology Project (SCP]_|22j and the High-z Supernova Search 
Team (HzST) and to the radio source data set of As far as the y 2 test is concerned, 
we shall show that the IRFP cosmology has a reduced x 2 value which is basically the same 
as that of the best-fit FRW solution. In the Appendix, we also apply the median statistics 



analysis 



2l| to the SNe la fit, and then use the Bayesian approach to test the robustness 



of the hypothesis that the correct model of the late Universe is the IRFP cosmology, and we 
show that the posterior probability that the late Universe is described by the IRFP cosmology 
is 36%. The corresponding probability for the best-fit FRW cosmology is only 24%. These 
numbers are obtained from a specific choice of the priors which was used in the literature 
already earlier, but of course, as always, it contains a certain degree of arbitrariness. 

The rest of this paper is organized as follows. In Section II we discuss the basic properties 
of the IRFP cosmology and introduce a possible extension of the model which describes the 
complete transition from a preceding FRW era to the IRFP cosmology. Section III is devoted 
to the calculation of the luminosity distance, and to the data analysis by means of the x 2 test. 
In Section IV we perform an analogous cosmological test, based on the angular diameter 
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distance-redshift relation. In Section V we discuss arbitrary power law cosmologies with 
variable G and A which generalize the IRFP model. In Section VI we derive a constraint on 
the IRFP from the CMBR acoustic peaks, and Section VII is devoted to the conclusions. 



II. THE MODEL 

A. The fixed point cosmology 

In this subsection we shall review the basic properties of the IRFP cosmology. It as- 
sumes spacetime to be homogeneous and isotropic so that it can be described by a standard 
Robert son- Walker metric containing the scale factor a(t) and the parameter K = 0, ±1 which 
distinguishes the three types of maximally symmetric 3-spaces of constant cosmological time 
t. The time evolution is governed by Einstein's equation R^ u — \g^ v R = —Ag^ + 8iiGT^ u 
with the conserved energy-momentum tensor T v = diag (—p,p,p,p). The equation of state 
is assumed to be of the form p(t) = wp(t), with w > —1 an arbitrary constant. Now we 
"RG improve" Einstein's equation by replacing G — > G(t) and A — ► A(t) which leads to the 



following coupled system of evolution equations 
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«V K 1 A 8tT _ 

- +- = -A + — Gp 3a 

a J a 2 3 3 

p + 3(l + w)-p = (3b) 
a 

A + 8np G = (3c) 

G(t) = G{k = k{t)), k{t) = A{k = k{t)) (3d) 



Eq. ()3a|) is the familiar Friedmann equation with a time dependent G and A, and eq. (J31 
expresses the conservation of T^. Eq. (}3*c|) is a novel integrability condition which ensures 
that the RHS of Einstein's equation has vanishing covariant divergence. 

The equations (j3d|) express the fact that the time dependence of G and A is the conse- 
quence of a more fundamental scale dependence of these quantities. We describe physics at 
a typical distance scale i = k^ 1 by means of a scale dependent gravitational action r^lg^] 
which should be thought of as a coarse grained free energy functional in the sense of Wilson. 
It encapsulates the effect of all metric fluctuations with momenta larger than the IR cutoff 
scale k, while those with smaller momenta are not yet "integrated out". In very general 
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terms, Tk describes the dynamics of fields "averaged" over spacetime volumes of linear ex- 
tension k~ x . The /c-dependence of Tk is governed by an exact RG equation whose precise 
nature is not important here. The quantum mechanical flow equation for the exact average 
action of QEG jj] would be one example, the classical equations of ref. are another. 
Generically Tj, is an arbitrary diffeomorphism-invariant functional depending on infinitely 
many generalized coupling constants which multiply all possible invariants which can be 
constructed from g^. We assume that we are in a regime where only the Einstein-Hilbert 
invariants J d A x^J —g R and J d A Xy/—g are important so that the only running couplings 
retained are their coefficients G(k) and A(k). If r& is of the Einstein-Hilbert form, the effec- 
tive, i.e. scale dependent field equations look like the conventional Einstein equation with 
the replacement G — > G(k), A — > A(k), where k \— > (G(k), A(k)) is an approximate solution 
to the RG equation. 

The final step of the improvement program consists in converting the /c-dependence of G 
and A to a time dependence. In ref. |l2| we discussed in detail that in a Robertson- Walker 
spacetime, to leading order, the IR cutoff should be identified as 

k(t) = | (4) 

where £ is a positive constant. 

Before discussing the solutions of the system (J3j) some definitions are convenient. We 
define the vacuum energy density the total energy density p tot and the critical energy 
density p cr i t according to 

. . A(t) , \ /r r \ 

PaW = 8lrG ( t y Aot(*) = P + Pa, (5) 

pUt) S 8^W) (l) (6) 

with H = a/a. Hence we may rewrite the improved Friedmann equation (jSaj) in the form 

a 2 + K 87T_,,. ,_ N 
— = -rG(t)p tot (7) 

CL o 

Frequently we refer the various energy densities to the critical density (J0J): 

Q M = — — , = , ^tot = + ^A = (8) 

Pcrit Pcrit Pcrit 

Another convenient abbreviation is fix = 1 — ^tot — 1 — — ^a- It follows from these 
definitions that 

A(t) = 3 Q A {t)H 2 {t) (9) 



and the Friedmann equation (JJJ) becomes 

K = a 2 [O tot - 1] = -a 2 n K (10) 

For a spatially flat universe (K = 0) we need p tot = p crit , as in standard cosmology. In this 
case the definition (jOJ entails 

Pu*(t) G(t) H 2 (t) = A (if = 0) (11) 

It is important to recall from Q| that the standard experimental value of Newton's 
constant, G Gxp , does not coincide with the value G(k = £/io) which is relevant for cosmology 
today, i.e. for t = to- In fact, G ex p is measured (today) at fc ex p — £/£, where i is a typical 
laboratory or solar system length scale. Thus, in terms of the running Newton constant, 
Gexp — G(k — fc ex p), since l/£ ^> l/t , and since in presence of several scales the relevant 
cutoff is always the larger one. On the basis of the antiscreening character of gravity we 
might expect that G(to) > G cxp therefore. 

Note that the critical density © is defined in terms of the cosmological Newton constant 
G(t). In particular in order to parametrize matter densities p obtained by non-cosmological 
measurements, it is sometimes more convenient to refer p to a modified "critical" density 
defined in terms of the experimental value G ex p: 

s S (12) 

A spatially flat universe requires p tot = p crit , but p tot might well be different from p e *f t . For 
the relative matter density referring to p e *f t we write 

V(t) = (13) 

Pcrit 

This definition implies that f2^ p = QuGe xp /G(t) or 

G(t) = TjSEpTTs G cxp (14) 

In [3] we investigated the consequences of the conjecture that, for k — > 0, the RG 
trajectory for the dimensionless Newton constant g(k) = k 2 G(k) and cosmological constant 
X(k) = A(k)/k 2 runs into an infrared attractive fixed point (<jr*, A*) with g* > and A* > 0. 
This means that the corresponding dimensionful quantities behave as G(k) = g*/k 2 , A(k) = 
\*k 2 for k — > 0. In the fixed point regime, the coupled system Q for K = and with 



k(t) = £/t was found to have an essentially unique attractor solution which governs the 
cosmology for t — ► oo: 

l/(3+3w) 



a(t) 



3^ 2 



t 4 /( 3 + 3 ^ (15a) 



P{t) = 9.(1+1)^^ (15b) 

G(t) = ^(l+ W )W (15c) 

o 

A(t) = 3|TT^F < 15d > 

Here the conserved quantity 

M = 8np(t)[a(t)f +3w = const (16) 

is the only free constant of integration. The system Q also fixes the parameter £: 

£ 2 = (17) 

? 3(1 + w) 2 K 1 ' 

The cosmology JT3J) yields p = p A = p CT it/2, p to t = Pent, or 

f2M = = 2' ^ tot = 1 (^) 

It is one of the most attractive consequences of the fixed point hypothesis that it leads 
unambiguously to these values of the relative densities. They are intriguingly close to, but 
not identical with, the values favored by the analyses of the recent observations within 
standard cosmology. We shall come back to this point later. 
The Hubble parameter of the cosmology (|15j) is 

H(t) = ~ (19) 
3 + 3w t 

and the deceleration parameter reads 

ad 3w — 1 , 

« s -f = — (20) 

We can use (II 9 j) in order to compute the age of the Universe, to, in terms of the present 
Hubble constant H(t ) = H : 

to = (21) 

° 3 + 3w K J 
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Clearly, for the late universe, the most plausible equation of state is p = 0, i.e. w = 0. 
In this case the fixed point solution describes an accelerated expansion a oc t 4//3 with the 
deceleration parameter q = — 1/4 =24|. For the age of the Universe we obtain t = f-f^T 1 - 
This is precisely twice the age one would obtain in standard cosmology with A = for the 
same value of Hq. 



B. Measuring G(to) and g*A* 

In the sequel we assume that the present Universe is described by the fixed point solution 
for K = 0, eqs. (|15p. and that this fixed point behavior started at a certain transition time 
t tr < to- Later on we shall determine t tr by postulating a plausible form of the RG trajectory 
running into the fixed point. 

In this section we discuss how, at least in principle, one can determine the present cos- 
mological Newton constant G(to) and the product g*\* which characterizes the fixed point. 
Note that the attractor solution (fTHjl depends only on this product but not on g* and A* 
separately. The experience with the ultraviolet fixed point of QEG 0, Q] suggests that the 
product <7*A* should be scheme independent (universal) while the factors g* and A* are not. 

Since the fixed point solution satisfies pit) = ptot/2, we can use ((TTj) in order to express 
G(to) in terms of observable quantities: 

16tt p(t ) 

By a very precise measurement of the Hubble constant and the matter energy density we 
can, in principle, determine G(to) and compare it to the laboratory value G e xp- It is clear, 
however, that before we can check whether G(t ) ^ G exp , the experimental situation must 
improve considerably. (Note also that the RHS of differs only by a factor of 2 from 
what one obtains in standard cosmology with A = 0, K = 0; in this case the prefactor is 

3/87T.) 

Because flyi(t) = 1/2 holds throughout the fixed point regime, we can interpret a devia- 
tion of G(to) from G ex p also in terms of a f2^ p (to) -value which differs from 1/2. From (|14jl 
we obtain 

G{to) = 2#) (23) 
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In order to get a first idea about the ratio G(to)/G exp we consider the hypothetical extreme 
case in which there exists no dark matter at all. Within the standard theory, the density of 
known forms of matter yields roughly 

Wo) > (24) 

Interpreting this figure within our model we find that G(to)/G exp ^ 50. Note that the ratio 
G(to)/G exp cannot be many orders of magnitude larger than unity precisely because the 
present density is so close to the critical one. 

As for the determination of g*A*, we start from the relation 

gA* = G(k)A(k) = G(t)A(t) = const (25) 

which, in the fixed point regime, holds for any cutoff identification k = k(t). Thus, g* A* = 
G(to)A(t Q ) where we may substitute 

A(*o) = \ Hi (26) 

which follows from (JUJ) with ^a(^o) — 1/2. Therefore we can express g*A* in terms of 
observable quantities either as 

g*K = \ G(t ) Hi (27) 

or as 

q /t4 

9m \ t = _L (28) 
y 32tt p(t ) y ' 

We see that we can determine (?*A* by measuring H and G(t ) or p(t ), respectively. Let 
us introduce the Hubble length £n(t) = l/H(t) and the Planck length £p\ = ^G exp , defined 
in the usual way in terms of the laboratory value of Newton's constant. In terms of these 
length scales (|2*7jl can be rewritten as 

3 G(t ) ( £ P} ' 2 



9 ' A - = 2^l4l)i (29) 

Since £ P i/£ H (t ) = O(10- 60 ) the product g*\* is of the order [G(t )/G cxp ] 10~ 120 . Most 
probably G(t ) and G exp do not differ by many orders of magnitude so that, roughly, g^A* = 
O(10- 120 ). 

We emphasize that in our scenario the smallness of this number does not pose any finetun- 
ing problem as it does in standard cosmology. In fact, g^X* is a fixed and universally defined 

10 



number which in principle can be computed from the RG equation. However, apart from 
being a difficult technical problem, this computation is possible only if we know the complete 
system of matter fields in the Universe. The number 10~ 120 reflects specific properties of 
this matter system coupled to gravity rather than an initial condition. 



C. A natural extension of the model 



Up to this point, our entire discussion was based upon the single hypothesis that, for 
k — > 0, the dimensionless couplings g(k) and X(k) are attracted towards a fixed point 
(<?*, A*) at g* > and A* > 0. By virtue of the system ® this information is sufficient in 
order to determine the cosmological evolution for t — > oo. In the K — - sector it was 
found to be given by the attractor solution (|TS|) . With this restricted knowledge about the 
RG trajectory, it is clearly impossible to determine how the Universe evolved towards the 
attractor, or to answer the question when the fixed point behavior set in. 

In order to have a model which applies also before the transition to the fixed point era 
we shall now generalize our hypothesis about the RG trajectory in a very "minimal", in 
fact the simplest possible way. We shall assume that the fixed point era of the cosmological 
evolution is preceded by an era during which G and A are approximately constant. 

To be precise, our first hypothesis is that the relevant RG trajectory, at least at a quali- 
tative level, can be approximated by 2 

, _ * Ik 2 for k < k tr 
G{k) = { a 1 (30) 
G< for k > k tv 

m = { x - k2 fOT k < ** (31) 

A< for k > k tT 

We fix the constant values G< and A< in terms of the "transition scale" k tr according to 



G< = (32) 
A< = \,kl (33) 



2 Eqs. (|3l)fl . could, for instance, be regarded as a model for the fi nal p art of a possible crossover from 
the UV to the IR fixed point, similar to what happens in 2D gravity 
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This guarantees that G(k) and A(k) are continuous (but not differentiable) at the transition 
point. Clearly a more realistic RG trajectory would be differentiable, and the transition from 
the regime with G, A = const to the fixed point regime would occur smoothly during a finite 
interval Ak tT centered about k tT . The trajectory (J3U|) . (|31j) should always be understood as 
an idealization of a smooth RG trajectory interpolating between G,A = const and the fixed 
point running. 

Our second hypothesis concerns the identification of the cutoff in terms of dynamical 
variables which in general could be of the form k = k(t, a(t), a(t), ■ ■ •). Guided by the 
discussion in we assume that the identification k = £/t, £ > 0, is valid for all k < k tT , 
where t tT = £//c tr denotes the time at which the transition occurs. (As we shall see, the scale 
k tr or the transition time t tT are not input parameters but rather are predicted by the model 
itself.) 

It is easy to solve the coupled system (jSJ) for t > t tr and t < t tr , respectively. For 
t > t tI , the most general solution is the one parameter family of attractor solutions fT5|) . 
with M. being the only free parameter. For t < t tr , the unique solution (with a(0) = 0) is 
the standard spatially flat FRW cosmology with a cosmological constant A< > 0. Its scale 



factor reads 



12 



a(t) = 

and the matter density is 



MG 



< 



2A< 



| cosh (1 + w) v / 3A7 t - lj 



-| l/(3+3w) 

(34) 



M 

P{t) = 8n[a(t)}^ (35) 
Note that the parameter M. has the same value in both branches of the solution because it 

is related to a constant of motion of the system Q by eq. (fTEJ) . 

As we discussed in detail in , the very existence of solutions to the improved evolution 

equations © is a nontrivial issue because for a given RG trajectory and cutoff identification 

this system is overdetermined. While in the case at hand it is very easy to solve the equations 

for t < t t r arid t > t tT separately, it is by no means guaranteed that, at t = t tv , the 

two solutions match in a physically acceptable way. Basically we are trying to adjust one 

parameter, t tr , in such a way that four functions, a(t) , p(t) , G(t) , A(i) become continuous at 

It is important here that in reality the transition from the classical FRW regime to the 
fixed point regime does not take place instantaneously, but during a finite interval of time, 
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At tr , which is centered about t tr . This transition period reflects the fact that the true RG 
trajectory is different iable, and that there is a smooth transition from G, A = const to 
G oc 1/k 2 , A oc k 2 during a finite but, by assumption, small interval Ak tr . The classical 
FRW solution is valid only for t < t tr — At tr /2, while the fixed point solution applies for 
t ^ t tr + At tr /2. During the transition regime 



the evolution is much more complicated, and we are not going to describe it explicitly. 
We assume that the system (jSJ) applies also in the transition period, and that a more 
realistic version of the RG trajectory (|30p. (j31j) gives rise to a continuous and differentiable 
solution for all four functions. The interpolating solution is likely to exist because during 
the transition the cutoff identification may well be much more complicated than k oc 1/t, 
the actual function k = k(t,a(t),a(t), ■ ■ •) being fixed to some extent by the requirement 
that (JHJ) is consistent. 

Thus our model provides an idealized description which ignores the transition regime, 
At tr 0. This is a sensible approximation if At tr is much smaller than t tr . Since in 
reality the functions a,p,G and A do have a certain variation between t tr — At tr /2 and 
t tr + At tr /2, it is clear, then, that in the model these functions should be allowed to have 
some moderate discontinuity at t — t tT - More precisely, we would consider it reasonable if 
a(t y t tr )/a(t \ ttr), say, is a small number close to unity. On the other hand, if this ratio 
should turn out many orders of magnitude smaller than unity, say, then this means that 
the Universe undergoes an enormous "inflation" during the transition regime. In this case 
a crucial qualitative element would be missing from our description if we just consider the 
classical FRW solution and the fixed point solution explicitly. Hence the idealized model 
would not be very useful probably. 

Let us now look at the properties of the extended IRFP model in detail. To start with, 
we assume that (|3*2|l and (|3l5|) hold true exactly so that G(t) and A(t) are continuous at the 
matching point. Eqs. and (jHHj) imply that 



This means that even if one allows for small discontinuities of G and A, the product G<A < 
is of the order of 10~ 120 . As the cosmological constant problem in its original form ('Why 



t tr - At te /2 < t < t t r + Attr/2 



(36) 



(37) 
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is A <§; mp^?" ) is related to the smallness of this number, we see that the mechanism which 
has achieved this tiny value must have been operative already before the FRW era with G < 
and A< has started. (Later on we shall assume that the FRW era began before the time of 
nucleosynthesis.) 

Using in = £/ k tT with £ given by (|17|) we obtain for the time of the transition: 



^ kl Z(l + w) 2 Kkl 3(l + w) 2 A< (38) 



Hence, in terms of A< or G 




(39) 



Another useful representation of t tT in terms of the age of the Universe can be obtained by 
eliminating g^A* from the second line of (J3~9"j) by virtue of (|27j) with H = 4/[3(l + w)t ]. 
The result is surprisingly simple: 



*» = W^y («) 

Now we are in a position to calculate a(itr) both from the fixed point solution a(t) = app(t) 
of ()15aj) and the classical FRW solution a(t) = apRw(^) i n (EH)- With (J3l?j) we obtain 

ofp(*«x)= — r-^ (4i) 
^< 

and 

".MG 1 1/(3+3u,) 

^FRw(^tr) = Cui — 7 = C w OFp(ttr) (42) 

where 

\ l/(3+3w) 

cosh(^)-l\ (43) 

In the following we continue the discussion for the most relevant equation of state p = 0, 
i.e. w = 0. 

In this case the ratio of the two scale factors, C w , has the numerical value Co ~ 1.55. It is 
encouraging to see that this number is indeed relatively close to unity. As Co > 1 the FRW 
scale factor at t = t tT is larger than the one of the fixed point regime. We discussed above 
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that, strictly speaking, the FRW solution cifrw(£) is valid only for t < t tr — At tT /2 and the 
fixed point solution app(i) applies only for t ^ t tT + At tT /2. During the transition a more 
complicated solution interpolates smoothly between czfrw an d ofp- The simplest possible 
behavior of the interpolating solution one could imagine is that a(t) stays approximately 
constant between t tr — At tT /2 and t tT + At tT /2. If one then chooses At tr such that 

a F Rw(Ur - Attr/2) = a FP (t tr + At tT /2) (44) 

the resulting scale factor is continuous and nondecreasing. Clearly, also a more complicated 
behavior of the exact solution is conceivable, but, even without knowing the details of the 
exact solution, the condition gives us a first idea about the size of At tT . Being interested 
in an approximate estimate only, it is actually more convenient to define At tr in a slightly 
different way, by aFRw(^tr) = &Fp(^tr + At tr ), which boils down to 

^ _ a F p(t tr + At tr ) _ ^ t tr + At tr ^j 4/3 ^ 

At tT _ ^3/ I 
Hr 

Even though At tr is not very much smaller than t tT itself, the ratio (|46|) indicates that there 
can very well be a window in which the idealized model is useful for a first orientation. 

Note that, as a consequence of the exact relation the discontinuity of pit) is neces- 
sarily of the same order of magnitude as that of a(t). 

From now on we shall be slightly more specific and assume that the FRW era with 
G = G < and A = A< has started before the time of primordial nucleosynthesis already 
(with w — 1/3 then). This assumption allows us to fix the value of G<. In fact, from the 
successful predictions of nucleosynthesis theory we know that the value of the cosmological 
Newton's constant during this epoch was very close to the laboratory value: G(t imc {) = G exp . 
Since, by assumption, G(t nuc i) = G < , this entails 

G< = G exp (47) 

With this identification, the transition time (j4()j) reads 

* = «M/|gy ( 48 ) 



OFp(^tr) \ t tr 

whence 

X/ Cq /4 - 1^ 0.39 (46) 
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This equation is quite remarkable. It relates the time when the fixed point behavior started 
to the ratio of the laboratory value of Newton's constant and the cosmological Newton's 
constant. Using (|23|) we can reexpress this ratio in terms of ^m P : 



t tI = t ^2QZ P (to) (49) 

Eqs. fi8|) and (}49|) are our main results of this section. 

It is convenient to characterize t tr by the corresponding redshift z tr . We have 

ago) ^oV /3 

and 

l + ^=(^) 2/3 = [2^o)]- 2/3 (51) 

Which values of z tT could we expect? The assumption f2^ p (t ) ~ 1/100 yields t tr > t /7 
and z tT ^ 12 , for instance. A distinguished period during the "recent" evolution of the 
Universe is the era of galaxy formation at about z ~ 3, say. If we speculate that the fixed 
point behavior started during this era we have 

z tI 3, t t r ~ (52) 

This scenario corresponds to f^ p (to) ~ 1/16. It has been shown in ref. that there are no 
growing "large" wavelength (which probably means super Hubble size) density perturbations 
in the fixed point regime. However, "small" -scale density perturbations can grow even in 
the fixed point regime. Hence z tr > 3 is not necessarily excluded by the theory of structure 
formation. 

We expect that if z tr is too small and the transition happened only recently, the fixed 
point solution cannot provide an accurate description of the present Universe because the 
transient phenomena of the transition region are still visible to some extent. Let us write 
At FP = t —t tT for the time which, according to the model, elapsed since the onset of the fixed 
point behavior. As a rough estimate we can say that the present Universe is already well 
within the fixed point regime, and the transient effects have become insignificant, provided 
Atpp ~ At tr - We rewrite this inequality as Atpp/^tr ~ At tr /itr- Then its LHS is given by 

^ = (1 + ^) 3/4 - 1, (53) 
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and 



yields Cq^ 4 — 1 for the RHS. Hence we are already well within the fixed point regime 



if 1 + z tr > C or 



z tT > 0.55 



(54) 



Thus we have a wide range of 2; tr -values where the model can be applied consistently. 

The estimate (J54j) is of a rather intriguing order of magnitude. In fact, the most distant 
supernova analyzed in 22] an d_j23j ] has a redshift of about z « 1, while the highest redshift 
bin for the radio sources in 28j] (see Section IIV|) is centered about z = 3.6. Thus, if the 
transition occurred relatively late, around z = 0.5, say, the current data might already tell 
us something about the transition from the FRW to the fixed point regime. This question 
will be discussed in detail in Sections IIHI and IIVI 



D. Summary of the extended IRFP model 

We investigated a simple cosmological model in which the IR fixed point regime of is 
preceded by a classical FRW era with constant values of G and A. While the hypothesis of 
this FRW era does not introduce any unknown new parameter, it allows for a determination 
of the time when the Universe entered the fixed point regime. In a nutshell, this time, 
t tr , obtains from the intersection point of the parabola Gpp oc t 2 with the straight line 
Gfrw(^) — G<- The transition time t tr turns out to be given by the ratio of the cosmological 
and laboratory value of Newton's constant or, equivalently, by the matter energy density 
^m P (^o)- If this density is of the order of the critical density, the ratio to/t tr is close to unity, 
i.e. the transition happened only "recently". In this sense the model solves the "coincidence 
problem" , the question why p and pa happen to be of the same order of magnitude precisely 
today. In our model we have p/p\ = 1 at any time later than t tT . 



III. TYPE IA SUPERNOVAE AND THE IRFP COSMOLOGY 



Recent supernova data from several sources (Supernova Cosmology Project |22j, High 
Redshift Supernova Team [23]) provide a valuable tool in cosmological hypothesis testing, 
since they contain information about the luminosity distance function di{z). In this section 
we shall apply the standard x 2 test in order to analyze the viability of the IRFP hypothesis. 
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In the Appendix, an alternative approach based on median statistics |2l| and a 
Bayesian model selection criterion is discussed in detail. 



A. The data set 

The data set used in all of our analyses consists of the 18 low redshift type la SNe from 
the Calan-Tololo survey 2(| plus 74 high redshift type la SNe, namely: 



42 SNe discovered by the SCP, w 
are published in Perlmutter et al. 



lose measured redshifts and apparent magnitudes 
221; 



• 32 SNe discovered by the HzST, who reported their measured redshifts and distance 
moduli in Riess et al. 

This gives a total of 92 supernovae. The oldest, most distant supernova is SN 1997ck, with 
a redshift of 0.97 = z m ^. 



B. Luminosity distance-redshift relations 



Let us assume that, at time t±, a distant galaxy emits light which is detected on Earth 
at to- The luminosity distance 4 to this galaxy is defined by the relation 

C 



T 



4nd 2 L 



(55) 



where £ is the absolute luminosity of the source and T is the measured flux. From the 
Robert son- Walker kinematics one obtains 

r to dt 
4 W). 

where 



a 2 (t 0/ 

d L = 77"^T" S 



a(t 



(56) 



sin x for K = + 1 
S(x) = { x for K = (57) 

sinh x for K — — 1 

Using a(to)/a(ti) = 1 + z we may rewrite the integral in (|5fi)) in terms of the redshift z: 



to 



dt 



fl 



dz 



, 1 



1 a(t) J H(z') 
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(58) 



(The subscript "0" on any quantity denotes its value at t = to.) In particular for K = 0, 



d L (z) 



l + z 



dz'- 



Hn 



H Jo H(z') 
In order to calculate H /H(z'), a further distinction is necessary: 



(59) 



1. The IRFP cosmology has K = and is given by (|15|) where we shall set w = from 
now on. The corresponding luminosity distance as a function of the redshift of the 
source reads 



drU) 



4(1 + z) 
Ho 



[l + z) l ' A -l 



(60) 



Eq. (p)Uj) is obtained either by inserting a oc 

£ 4/3 

into (JSTIj) and eliminating t± in favor 

of z, or by using (J53J) with 

tf(£)=tfo(l + ^ 3/4 (61) 

which follows from (fTTJj) . 



2. In standard cosmology one has the familiar relationship 



which can be used to derive 
(1 + *) 



1 + zfQuo + Ovo + (1 + zfVL 



K0 



1/2 



(62) 



driz) 



Hoy/P 



K0\ 



A"0 



i + z') 3 n M0 + n A0 + (i + z'fn 



AT) 



-1/2 



(63) 

As for the actual form of the function S, eq. (fTU|) shows that the cases K = +1, and 
—1 are realized if fixo < 0, ^ko — and f2 K o > 0, respectively. 

3. In the extended IRFP model introduced in subsection III CI we combine a K = FRW 
cosmology, valid for t < t tr (z > z tT ), with the fixed point cosmology which applies for 
t > t tT [z < z tT ). For arguments z < z tT , the corresponding function d^{z) is given by 
(|50|) . For arguments z > z tr we divide the z'-integral of in a IRFP plus a FRW 
piece: 



driz) 



l + z 
Ho 



dz' 



Ho 



H(z> 



+ I dz' 

IRFP 



Ho 



H(z> 



FRW 



(64) 



The integrand of the IRFP piece is given by (jBTjl and for the FRW term we use (16 2 j) 
with fi K o — 0: 



driz) 



l + z 

Hn 



ZtT 



dz'(l + z')- 3 / 4 + / dz' (1 + z') 3 Q M0 + n 



'AO 



-1/2 



(65) 
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The quantities f2 M0 and ^ao = 1 — ^mo are n °t the physical densities at t — to, these 
are I^mo — ^ao — 1/2? but rather fictitious values which would result from using the 
standard FRW differential equations also for t between t tv and to- The quantities Quo 
and f^Ao are determined by the requirement that the FRW- and the IRFP-solution 
match at z = z tT ; they depend on z tr therefore. 

In the present context it is most convenient to define t tv by the requirement that 
Qu{t) and are continuous at the transition. 3 Then the specific FRW theory 

which connects to the fixed point cosmology is singled out by the requirement 

nT w (t tr ) = nf FP (t tT ) = 1/2 



and similarly for Qu- On the FRW side we can use the general result 4 

-l 



AO 



'i + zYn 



MO 



n 



AO 



in order to relate f2^ RW (t tr ) to Q\ . Letting z 
Quo — > ^mo, eq. yields 



1 + z) 2 Q K0 



Ztr, ^AOtr) = 1/2, Q 



AO 



^ao — (1 + z tr ) 3 Q M0 



which, with Q 



MO 



1 — f^Ao, leads to 



n 



MO 



1 + (1 + Ztr 



q 



Ztr 



AO 



(66) 
(67) 

(68) 

(69) 
(70) 



1 + (1 + ^ tr ) 3 

Using (J69|) and (|70|) in (J65)) we obtain the final result for the luminosity distance in 
the case z > z tT : 

l + z 



d L {z) = i±i{4(l + ^ r ) 1 / 4 -4 



+ 



1 + (1 + Ztr 



1/2 



dz 



[l + z 'f + ll + Ztr f 



-1/2 



(71) 



It is customary to express luminosity distances in terms of the distance modulus ^o 
m(z) — M according to 

dAz 



m(z) = M + 5 log 10 



1 Mpc 



25 



(72) 



3 The ttr thus defined coincides with the one from subsection III CI within Ai tr . 

4 Eq. jHZt follows from by noting that p\ = const in the standard case, and that p C rit oc H 2 . 
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where m and M are the apparent and absolute magnitude, respectively. Measuring dt in 
units of the present Hubble radius, l/H , the dependence of the apparent magnitude on the 
"Hubble-free" luminosity distance H cIl reads 

m( z ) =M + 5 log 10 (h d L {z)) (73) 

where 5 

M = M - 5 log 10 (H ■ 1 Mpc) + 25 (74) 

is constant within an ensemble of standard candles. 

The HzST results are presented in terms of the distance modulus, while the SCP published 
the estimated effective i?-band magnitude which relates to the HzST results through 

mf = M B + no (75) 

where Mb is the peak 5-band absolute magnitude of a standard type la supernova. We 
chose to translate the HzST data from distance moduli into apparent magnitudes, leaving 
the SCP data unchanged. 

In the expression for the apparent magnitude, Mb and Hq appear only in the combination 
M. B = M B — 5 log 10 (-Ho ' IMpc) + 25. We have used the value of Mb that is obtained from 
the analysis of the low-redshift SN sample of the Calan-Tololo survey 
the expression (|75|l for the apparent magnitude becomes 



alone. If z <C 1, 



m B (z) = M B + 5\og w z (76) 

Fitting the measured apparent magnitudes to this formula yields .Mb = —3.32. 

The total data set now consists of 92 measured apparent magnitudes m B ieas , henceforth 
denoted m™ eas , i = 1, • • • ,92, which must be compared to the theoretical expectation for 
the Hubble-free expression for the apparent magnitude, 

m B (z) =M B + 5 log 10 (# d L (zf) (77) 



Fig. displays the apparent magnitudes m™ eas as a function of the corresponding mea- 
sured redshifts z { . The error bars on the data points indicate the statistical errors (jj of 

5 There should arise no confusion with the invariant with the same name M. which is also customary. 
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FIG. 1: The measured apparent magnitudes of the supernovae as a function of their redshift. The 
continuous line represents the prediction of the IRFP cosmology, the dashed one is the best-fit 
FRW model, and the dot-dashed line is a flat FRW model with zero cosmological constant. 

m meas estimated by SCP and HzST, respectively. Fig. ^ also shows the theoretical curve 
f!77|) with the function d^z) pertaining to three different cosmological models: the IRFP 
model (without a FRW part), a FRW model with Q M0 = 0.3, ^ao — 0.7, and a FRW model 
with £7mo — 1, ^ao — 0. The second model is the one which, in the framework of standard 
FRW cosmologies, fits the data best. The question is whether the IRFP leads to an even 
better fit. 



C. x 2 analysis of the currently available 

The standard x 2 test is the most common 
between the predictions of some cosmological 



data 

choice in order to quantify the agreement 
model and the supernova data. From the 
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triples (z i ,mf ieas ,ai), i = 1, • • • ,92, provided by the dataset described in IIII Al we have to 
compute 



where m* heor = m^,{zi) is the theoretical value predicted by eq. (|77|) with the function 
(Il(z) computed within a specific model. As for the fixed point model, we must distinguish 
the cases z tT > z max and z tr < z mSuX . In the former, the Universe had already entered the 
fixed point regime by the time the oldest supernova exploded; in the latter, some of the 
supernova events took place during the preceding FRW era, and their photons detected by 
our telescopes have witnessed the transition to the fixed point regime. 

If z tT > z max then di,{z,i) is given by (plj) for all supernovae. The evaluation of the sum 
(|T8j) yields a value of y 2 per degree of freedom (x 2 /dof) of 1.59. In Table |U we compare this 
number to the corresponding \ 2 values for three competing FRW cosmologies: 

• The best-fit FRW model with vanishing cosmological constant and Quo anywhere in 
the interval [0, 3]; 

• The best-fit FRW model with Qmo £ [0,3] and ^ao G [—1,3] such that Qmo + ^ao = 1 ; 
i.e. spacetime is spatially flat; 

• The best-fit FRW model with Q M0 e [0, 3] and f^o £ [ — 1, 3] arbitrary. 

We see that the x 2 /dof value of the IRFP cosmology is about the same as those for 
the FRW models. Taking the numbers in Table |J at face value, the fixed point model 
performs even slightly better than its FRW competitors. This is a very encouraging result, 
in particular since, contrary to the FRW models, the fixed point model has no free parameters 
which could be adjusted. 

Up to here we assumed that z tr > z max so that even the oldest supernova in the dataset 
would not "know" about the onset of the fixed point epoch. Assuming the opposite case 
Ztr < Zm&x now, some of the older supernovae exploded already during the FRW era which, 
according to the extended IRFP model, preceded the fixed point epoch. The function di(z) 
of the extended IRFP model depends parametrically on z tr . For supernovae with < z tr 
the value of d L (zi) is given by equation (|5Uj) . while for Z\ > z tr eq. (fTTj) must be used. As 




(78) 
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a consequence, y 2 is a function of z tr now, and in principle one could hope to extract the 
value of z tr from the data by looking where x 2 = X 2 { z tr) is minimum. In the ideal case the 
function x 2 { z tr) would have a pronounced minimum at some value of z tT , and we would then 
be entitled to conclude that the fixed point regime started at this redshift most probably. 

Using again the data set of subsection 1111 Al we obtain the function \ 2 — X 2 { z tr) which is 
plotted in Fig. |21 Obviously there is no strongly preferred value of z tr - There seems to be a 
minimum at z% t ~ 0.18, but it is too shallow to derive any statistically significant conclusion 
from it. 

Thus we must conclude that the statistical quality of the currently available supernova 
data is not good enough in order to discriminate between the fixed point model proper 
(without a FRW era) and the one-parameter family, parametrized by z tr , of extended fixed 
point models. 

D. x 2 analysis of Monte Carlo data for SNAP 

It is an important question whether it would be possible to detect a "smoking gun" of the 
transition to the IRFP regime in the supernova data if a better statistics were available. In 
fact, the SNAP (Supernova Acceleration Probe) satellite is aimed to improve this statistics 



TABLE I: The value of x 2 P er degree of freedom for several cosmological models (the errors on the 
best-fit parameter values correspond to the 1-a confidence region). 



Model a 


X 2 /dof 


Parameter values 
at minimum 


Fixed point model (z tr > z max ) 


1.59 




Extended fixed point model 


1.60 


z . _ n 17 +0.80 
«tr — U-^'-O.lO 


Best-fit FRW with fl A0 = 


1.62 


n M0 = 0.00+ ' 05 


Best-fit FRW with ft A o + ^mo = 1 


1.60 


Quo = 0.40±°;J° 


Best-fit FRW 


1.62 


Quo = 0.25l°;|j, n AO = OAot°oll 



a The variation intervals for CImo and Qao i n the fits are [0,3] and [-1,3], respectively. The transition redshift 
goes from to 0.97, the redshift of the farthest SN in the set. 
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FIG. 2: The x 2 value of the extended IRFP model as a function of the redshift ztr obtained from 
the currently available data set. 

at a significant level. Its goal is to obtain a super data set more than one order of magnitude 
larger than the currently available one, with an improved control over systematic errors, to 
redshifts up to about z ~ 1.7. We therefore decided to simulate the expected results in 
order to forecast the impact of this enlarged data set on the constraints for z tr . In particular 
we used the SNOC code, developed by A.Goobar et al. j^J, which takes into account 
gravitational interactions (lensing) and extinction by dust, both in the host galaxy and 
along the line-of-sight, and modified it by implementing our new distance- redshift formula. 

The modified SNOC code based upon the extended IRFP model was used to generate 4 
different Monte Carlo data sets, each consisting of 2000 type la supernovae in the redshift 
interval [0.1, 1.8]. The 4 simulations differed with respect to the transition redshift z^ mul of 
the underlying IRFP model; we chose z^ mul = 0.2, 0.4, 0.6 and 1.0, respectively. We then 
performed a x 2 - ana lysi s of the 4 data sets as we did with the "real" data in the previous 
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FIG. 3: The function x 2 { z tr) from a simulated data set of 2000 supernovae for various values of 
the transition redshift z^ mul . 

subsection. Fig. Elshows the resulting x 2 as a function of z tT (not to be confused with ). 
We find that, if the transition took place late enough, at a redshift of 0.2, say, there is a clear 
minimum in the x 2 { z tr) curve, and it should be possible to determine z tT from the SNAP 
data 6 . On the other hand, if it occurred much earlier, at a redshift of 1.0, say, the x 2 { z tr) 
curve does not prefer any specific value of z tr . For z tT > 0.5 this curve is almost flat. This is 
related to the fact that, if z tT ^ 1, the luminosity distance functions d^z) of (|6TJ|) and (|7T|) 
are virtually identical for arguments z G [0,2], say. 

The lesson we learn from this Monte Carlo investigation is that, if we improve the statistics 
by a factor of about 10, there is a significant hope of either actually detecting the transition 
to the fixed point regime if z tr ~ 0.5, or at least of putting a lower bound on z tr if z tr > 0.5. 

In the Appendix and in Section V we resume the analysis of the currently available 

6 In view of the discussion which led to the estimate l|54(l it is clear, however, that a quantitatively correct 
description of this late transition would require a more sophisticated model than that of Section [n] 
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"real" supernova data. Since they are anyhow not sufficient to determine z tI we assume that 
zti > z max there. 



IV. COMPACT RADIO SOURCES AND THE IRFP COSMOLOGY 

In the previous section we have shown how a class of astronomical objects of well-known 
intrinsic properties can help probing the cosmic evolution via the luminosity distance-redshift 
relation. 

A similar, yet independent test is based on the angular diameter distance-redshift relation. 
The angular diameter distance dA{z) of an object with redshift z and proper diameter / which 
subtends an angle 6 as seen by a terrestrial observer is defined as 

dA{z) = = (79) 

where we have introduced the characteristic angular size D = IH which is to be interpreted 
as an angle given in milliarcseconds (mas). By measuring the object's redshift and #-angle 
one can compute its value of d^{z) and compare it to some particular cosmological model, 
provided its diameter I is known. By comparing the predictions of several different models, 
the one which best fits the data can then be estabilished. Of course, a class of standard rods 
(objects with the same characteristic extension I) is needed in order for this determination 



of dA{z) to be viable. Recently 28 1. with some manipulation and redshift-binning, a set of 



330 compact radio sources has been applied to this purpose, thus giving the opportunity to 



constrain the free parameters o: 
cardassian expansion scenarios 



several cosmologies, from FRW and power-law models to 
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This section intends to address precisely the issue of confronting the IRFP model with 
the new data. The statistical framework is analogous to the one adopted in 11111 both the 
proper and extended versions of the IRFP model are compared to FRW cosmologies via the 
X 2 test. The x 2 t es t for the extended IRFP model is applied to the data set in order to give 
an estimate of the transition time t tT . 

A. The data set 

The class of standard objects to be used as distance indicators is compact (< 100 pc) 
radio sources. Specifically, the set consists of 330 sources with redshifts Zi in the interval 
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[0.011,4.72] and angular diameters 9i, extracted from 5 GHz VLBI contour maps in the 



literature, as reported by Gurvits et al. 
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Several reasons point out why these are likely to be standard objects: first, evolutionary 
effects are expected to be negligible with respect to galaxies or extended two-lobes radio 
sources, since the characteristic time scale of the former (some ten years) is very tiny com- 
pared to the age of the universe, or to the Hubble time H^ 1 . Second, while the physics 
of extended systems is more likely to be influenced by the different properties of the in- 
tergalactic medium encountered at different redshifts, the features of compact sources can 
in principle be described by a few parameters, like the mass of the central black hole, the 
magnitude of the magnetic field and of the angular momentum, and a few more. 

Three tasks have thus to be completed before the data set can be applied to the cosmo- 
logical analyses: 

(i) To control and put limits on the few basic parameters governing the radio source, so 
as to have a collection of objects which is as uniform as possible; 

(ii) To homogenize the set, which has been assembled from contour maps published by 
different collaborations, rather than from a unique set of observations; 

(iii) To investigate possible evolutionary effects (linear size - redshift dependence), as well 
as linear size - luminosity and linear size - spectral index dependences, and minimize 
their influence on the cosmological implications. 

To these aims, Gurvits at al. restrict the set by retaining only those sources with 
luminosity Lh 2 > 10 26 W/Hz and spectral index —0.38 < a < 0.18. Since a moderate 
dependence of a on the angular size seems to be estabilished (see Fig. 7 in the last 

relation has the further effect of constraining the dimensions of the sources. 

The new group of 145 sources is then binned into 12 redshift intervals, and a median 
redshift and angular size is calculated for each one, together with the standard deviation 
<7j. This procedure has the additional effect of reducing the influence of particularly large 
objects in the set, and providing simple means to deal with unresolved sources (see cases J, 
L and S in 2^| for details). 

The final data set consists of 12 triples (zi,8i,ai) in the redshift interval [0.52,3.6]. The 
data points, together with the three curves representing the proper IRFP prediction and the 
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FIG. 4: The 12 data points (with corresponding uncertainties) and the theoretical predictions of the 
proper IRFP model and of the FRW models with = 0.3, = 0.7 and with Om = 1, = 0. 
The angular size D is assumed to equal its best-fit value D = 1.30 . 

FRW models with Qu = 0.3, Q\ = 0.7 and with Qm = 1, = 0, are presented in Fig. |U 



B. The angular diameter distance-redshift relation 



From (J79|) . following simple reasoning, one can easily show that in any spatially flat FRW 
spacetime 



d A (z) 



driz) 



[1 + zY H (l + z) 



dz' 



i H 
H(z> 



(80) 



In IIII Bl we demonstrated that, in the extended IRFP model, for z greater than z tr , 



dz'- H ° 



/o H(z') 
while, for z < z tT 



4(l + z tr ) 1 / 4 -4+ l + (l + z tr ) 



1/2 



dz' 



-1/2, 



dz'—r-T 
o H(z') 



(l + z) 1 / 4 -! 



l + z'f +{l + Z tx ) 6 (81) 



(82) 
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It is evident from the above formulas that (Ia(z) also depends on the transition time, as long 
as the transition does not occur at a redshift zt T which is greater than the maximum z of the 
sources in the data set (,2 max ~ 3.6, see below.). We can therefore investigate the possibility 
of a transition to the fixed point regime in the redshift interval [0, 3.6]. This window is much 
larger than the one allowed by the supernova set adopted in 11111 



C. x 2 analysis and determination of the transition time 

A x 2 function can thus be constructed in the usual manner: 



(83) 



»=i 

Here 0j stands for the observed values of the angular size with errors <7j, and 

/ D 



d A (zi) H d A (zi) 

is the theoretical value. Let us note that the x 2 value is also a function of the "Hubble free" 
diameter D of the source (it does not depend on H though; see Eq. (j80|) for (1a(z)). The 
neatest way to solve this problem is to calculate a 2-dimensional x 2 i n the parameter space 
of Ztv and D, and then marginalize over D to make the results independent of its value j^. 
The level contours for the probability P(z tr ,D), associated to x 2 { z tr,D) are shown in Fig. 

P(z tT ,D) is maximum at (z tr = 0.08, D = 1.3 mas), where x 2 has a corresponding value 
of 1.86. In Table ITT1 this result is compared to other cosmologies. Marginalizing over D, i.e. 
defining a new probability as 

~ s f 2 dDP(z tI1 D) , 
max [P(z tr , D)\ 

one obtains the probability distribution for z tT alone, as shown in Fig. |U1 

Clearly, the situation is even more problematic than for SNe la: the probability P(zt T ) 
hardly seems to vary, and a determination of z tT is by no means possible with the current 
data. However, again we may conclude that the IRFP model proper performs as well as 
standard cosmology. 
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V. BEYOND THE IRFP MODEL: GENERAL POWER LAW SOLUTIONS 



In this section we are going to put the fixed point cosmology into a more general perspec- 
tive by "embedding" it into a broader class of cosmologies with a time dependent G and A. 
Let us give up for a moment the idea that the t-dependence of G(t) and A(t) stems from 
an underlying RG trajectory via some identification k = k(t). Then we are left with the 
system of equations ()3al3bl3c|) without pH]) . This system has been studied in the literature 
already long ago j^lj]. The problem is that (j3a!3bl3c|) is underdetermined, so that in order 
to obtain a unique solution one has to "invent" some additional condition on a, p, G and A 
on an ad hoc basis. For instance, without providing a deeper explanation, it was assumed 
3l| that Newton's constant follows a power law G(t) oc t n with an arbitrary, not necessarily 
integer exponent n. With this additional equation, the system (j3al3b!3c|) has the following 

TABLE II: The value of x 2 per degree of freedom for several cosmological models (the errors on 
the best-fit parameter values correspond to the 1-a confidence region). 



Model a 


X 2 /dof 


Parameter values 
at minimum 


Fixed point model (zt r > z max ) 
Extended fixed point model 
Best-fit FRW with Q A0 = 
Best-fit FRW with ft A o + ^mo = 1 
Best-fit FRW 


1.88 
1.86 
1.91 
1.86 
1.91 


n _ i 47 +o.06 
~ _ n nR+ - 80 n - 1 qn+ - 30 
O M0 = 0.55+° 20 , = 1.20±g;i8 
O M0 = 0.60j;^, J D = 1.20l°|° 
^mo = 0.151^1, n A = 1.201^, D = 1.8G±oj° 



The variation intervals for Omo and £Iao in the fits are [0,3] and [-1,3], respectively. The transition redshift 
z tr goes from to 4. 
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FIG. 5: Contour levels for the probability P(zt r ,D). 



solution for K = 0: 



a(t) -- 

P(t)~- 
G(t) 
A(t) 



3(1 + w) 2 
2(n + 2) 

(n + 2 



M C 



1/(3+3to) 



t 



(n+2)/(3+3to) 



12tt (1 + w) 2 C t n + 2 
Ct n 

n(n + 2) 1 



(85a) 

(85b) 
(85c) 
(85d) 



3(1 + w) 2 t 2 

Actually ()85)1 describes a 2-parameter family of solutions labeled by the parameters Ai and 
C. The fixed point cosmology (fTB]) is the special case of which is obtained by setting 
n = 2 and C = 3(1 + u>) 2 g*A*/8. In fact, the exponent n = 2 is an unambiguous prediction 
of the fixed point hypothesis. It obtains not only for the simple cutoff identification k oc 1/t 
but even for an arbitrary function k = k(t). This can be seen as follows. On the one hand 
we have, in the fixed point regime, G(t)A(t) = G(k)A(k) = g*\*; on the other hand, eqs. 
x| and ()85d|) yield G(t)A(t) = const ^ if, and only if, n = 2. 
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FIG. 6: The marginalized probability P(z tT ). A shallow maximum can be seen at zt r = 0.16. 

The luminosity distance for the models (|K5|) is easily worked out. It is independent of Ai 
and C, but it does depend on the exponent n. For a power law expansion a(t) oc t a we have 
in general 

- (^t) ^ K 1 + ~~> , - 1/ " - !] < 8 «) 

and (|85ajl yields in particular 



a = (87 

3 + 3w K J 

We used the currently available supernova data in order to test the cosmologies (|85j) . 
Proceeding as in Section IIII C\ we performed a x 2 analysis where 77j,* heor was computed from 
the luminosity distance (jSHJ), (jHTjl with w = 0, i.e. for a = (n + 2)/3. Fig. [3 displays the 
resulting y 2 values as a function of the exponent n. Quite remarkably, we observe a clear 
minimum of \ 2 at n = 2, which is precisely the exponent predicted by the fixed point model. 
Stated differently, in the space of all cosmological models with a power law expansion a oc t a 
it seems to be exactly the t 4//3 -expansion of the IRFP model which fits the supernova data 
best. 
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FIG. 7: The x 2 values of the general power law models as a function of the exponent n. 
VI. CONSTRAINTS ON z tr FROM THE CMBR FIRST ACOUSTIC PEAK 

The recent data from WMAP have determined the position of the first acoustic peak in 
the CMBR spectrum with great accuracy |22|. The spectrum depends on the complicated 
interplay between spacetime geometry and microphysics of the recombination era but the 
position of the first Doppler peak is essentially a measure of the angular size of the sound 
horizon, A8 S , at the recombination epoch (z rec « 1300). One has 

1 



d A (z Tec ) (gg) 



dfj(z rec 



A6 S 

where dA(z rec ) is the angular diameter distance to the surface of last scattering, and dfj(z Tec ) 
is the proper radius of the sound horizon at decoupling j^. The precise proportionality 
factor on the LHS of (J%%|) is dermined by the microscopic model of the recombination. In 

rom ~ 0.85° for a simple A = model, to about 
. Thus, strictly speaking the A^ s -value inferred 
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standard FRW cosmology A9 S ranges 
~ 0.60°, in a more refined calculation 



from a given measured multipole number i\ is model dependent, but the model dependence 
changes A6 S at most by a factor of order unity. Taking it for granted that A8 S is about one 
degree or slightly less, i.e. of the order of magnitude predicted by standard cosmology, we 
can derive an upper bound on z tT in the following way. 

The key observation is that both dA{z rec ) and d s H (z vec ) tend to increase when we replace 
a piece of FRW evolution by the FP evolution. For z tr small enough, their ratio A6 S stays 
almost constant. 

As for dA(zrec), Eq.(80) with (pjj) provides us with (1a(z) for a power law expansion 
a(t) oc t a . We can use this formula for a comparison of the t 4//3 -FP-evolution and the t 2//3 - 
FRW-evolution where we neglect the cosmological constant for a first orientation. In the 
extreme case of a transition during the recombination 6r£l ; Zfa — ^rec 

= 1300, one finds 

d A (z 



« 10.3 d A (z Iec ) (89) 

a=4/3 a=2/3 

We conclude that for realistic redshifts z tr <C z TCC the value of dA(z TCC ) is bounded above by 
about 10 times its value in standard cosmology. 

As for djj(z iec ), we recall that in any (Qmo, ^ao)-FRW model the sound horizon at redshift 
z ^> ^MO) ^ao has proper radius 

4(*)« 77^(1 + *)~ 3/2 (90) 
"Ov " 'mo 

where c e ff is an effective speed of sound. In the extended FP model, the FRW cosmology 
preceding the FP era is characterized by the densities f^Mo and f^Ao of Eqs.(69), (70). This 
allows us to compare dfj(z rec ) in the extended FP model to the corresponding value in the 
Q uo = 1, a = 2/3 FRW model: 



dff(z TeC/ 



= + + ^tr) 3 4(*rec) (91) 

FP o=2/3 

Combining (|55jl to (|9~T|) we observe that the increase of dA(z rec ) certainly cannot compensate 
for the increase of d s H (z mc ) so as to keep A9 S unaltered if the square root in (pTTj) is larger 
than 10.3. The condition yjl + (1 + Ztr) 3 < 10.3 leads to un upper bound for the transition 
redshift: 

ztr < 3.7 (92) 

This is a remarkable order of magnitude. It indicates that the transition should have taken 
place at or after the era of structure formation. This conclusion is in accord with the picture 
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arising from the analysis of density perturbations ^| . The actual value of is much smaller 
probably. For a more precise estimate the details of the microphysics should be taken into 
account properly which is beyond the scope of the present paper. 

VII. CONCLUSIONS 

In this paper we used the SCP and HzST data on high redshift type la supernovae and 
the radio source data from Gurvits et al. 1281 in order to test the infrared fixed point model 

n 

proposed recently 16J. It predicts a time dependent cosmological and Newton constant 
whose dynamics arises from an underlying renormalization group flow, and it leads to a 
characteristic t 4//3 -time dependence of the scale factor in the late Universe. The latter gives 
rise to a very particular distance-redshift relation which, in principle, can be verified or 
falsified using distant supernovae as standard candles and compact radio sources as standard 
rods. 

The results of our investigation based upon the \ 2 test can be summarized by saying 
that the IRFP model performs at least as well as the best-fit FRW model in reproducing 
the supernova and radio source data. Moreover, contrary to FRW cosmologies, it has no 
free parameters and thus actually explains why flu an d are of the same order today. In 
fact, in the flat case, they are predicted to be exactly equal (Qm = = 0.5), today and at 
any time in the future. 

In the space of all cosmologies with a power law expansion a oc t a , it seems to be precisely 
the exponent a = 4/3 predicted by the IRFP model which yields the best fit to the data. 

In this paper we also extended the original IRFP model of 16] by matching the fixed point 
regime with a preceeding FRW era. We found that the statistical quality of the currently 
available data is not sufficient in order to determine when the Universe entered the fixed 
point epoch. However, using a Monte Carlo simulation we saw that the data expected from 
the forthcoming SNAP mission can teach us something about the transition to the fixed 
point regime provided it happened late enough, at a redshift smaller than 0.5, say. 

In conclusion we think that the fixed point model is a viable alternative to the quintessence 
scenarios and certainly deserves being studied further. In particular it would be interesting 
to compare in more detail its predictions for the microwave background radiation to the 
data. We hope to address this issue in a future publication. 
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APPENDIX A 

In this appendix we resume the analysis of the supernova data in terms of the IRFP 
model proper, using a different statistical method: median statistics and a Bayesian model 
selection criterion. 

1. Median statistics 

Median statistics j3, I2]] is more easily implemented than the x 2 analysis and also less 
demanding with respect to the assumptions about the data because uncertainties need not 
be normally distributed, and no prior knowledge of measurement errors is required. 

According to median statistics, the operation of calculating the likelihood of a particu- 
lar model boils down to counting: one simply enumerates how many supernovae are, say, 
brighter than expected. Now, the probability that, when no systematic errors are present, 
n out of 92 SNe are brighter than expected (92 — n being, of course, fainter) is given by 
the standard binomial distribution P(n, 92). As for the fixed point model proper (without 
a FRW epoch) and the data set of IIII Al this number is found to be to n = 53 which has 
binomial probability P(53, 92) = 0.029. For comparison, we report in Table ITTT1 the same 
analysis for other models. Notice that, if a model contains one or more free parameters, in 
column 1 we refer to those regions in the parameter space where the binomial probability is 
maximum. 

At first sight the IRFP model seems to perform rather poorly in comparison with the 
FRW models. However, it is important to understand that the figures in the last column 
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TABLE III: Binomial likelihood after 92 SNe. (For models with free parameters, the numbers 
pertain to those regions in parameter space where the probability is maximum.) 

Model Brighter Fainter Binomial likelihood 



Fixed point model (z tr > -2 max ) 53 39 0.03 

Extended fixed point model (z tl < 0.01) 46 46 0.08 

FRW with A o = 0, < n M o < 0.1 40 52 0.04 

FRW with fi A0 + n M o = 1, 0.4 < n M o < 0.5 46 46 0.08 

FRW with fi M0 and ft A0 in the shaded region of Fig. El 46 46 0.08 



of Table ITLTl are not a measure for the correctness of the various models. In fact, the FRW 
models contain free parameters which always can be adjusted such that there is about the 
same number of brighter and fainter supernovae (at least if one allows for a cosmological 
constant). This is not possible with the IRFP model which contains no free parameters. 
Hence the naive median statistics is necessarily "unfair" towards the fixed point model. The 
correct way of assessing the "relative correctness" of the various models with their different 
number of free parameters is provided by Bayesian statistics, to which we shall turn in the 
next section. 

2. Bayesian model selection 

a. The method 

In the Bayesian model selection the data set DS is assumed to have arisen from one 
of several possible models (or hypotheses) Mi, ■ ■ ■ ,Mn- An a priori probability -P(Mj) is 
assigned to each model in order to measure the likelihood of Mi when no other information 
(e.g., observational) is available. This "prior" depends exclusively on the way each model is 
structured: we shall see below how to relate this number to the number of free parameters 
it presents. Let P{DS\Mi) be the likelihood that a data set DS, such as the one that 
is actually observed, is attributable to Mj. Then, following Bayes' theorem, the posterior 
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probability P(Mi\DS) that the model Mj is responsible for the observed data set DS is 

P(MWS) - P{DSm Pm (Al) 
P{M % \Db) - ^ p{DS \ Mj) p{Mj) ( A1 ) 

The likelihood P(DS\Mi) equals the statistical probability that, given Mi is true, a set of 
n observations would result in the measured outcome DS. Therefore, these numbers are 
just the experimentally determined binomial probabilities that we have already discussed in 
section IA 11 



b. Assigning a priori probabilities 

Let us now assume we are trying to determine the "degree of belief" P(Mi) that model 
Mi is the correct one, before we know anything about the way it reproduces observations. 
Apart from the constraint that 

N 

£P(M 4 ) = 1 (A2) 

i=i 

we have no other clues. Basically, we would like to quantify the extent to which a given 
model is constructed on an ad hoc basis. According to Ockham's razor principle, the less free 
parameters are present, the more a model can be regarded as physically plausible. Bayesian 
priors can effectively help implementing this principle if one relates the initial probabilities 
P(Mi) of a model to the number of free parameters it contains. Following j^, we then set 

where N$ is the number of free parameters of M$. It is easy to see that this definition 
guarantees (|A2|) . provided we share the corresponding probability into equal parts whenever 
there is more than one model with the same iVj. 

In this paper we consider the following hypotheses: 

1. Mi: The IRFP cosmology correctly describes the cosmic evolution for z < z tr , with 

^tr ^ 2 maX) 

2. M 2 : The Universe is described by a zero-A FRW cosmology, with f^Mo > 0.1; 

3. M 3 : The Universe is described by a flat FRW model, with 0.4 < Q M0 < 0.5; 
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4. M^. The Universe is described by a general FRW model, with £7mo and £7ao both in 
the shaded confidence region indicated in Fig. |HJ 

The number of free parameters and the corresponding priors are listed in Table IIVI (Here 
and in the following "region 1" stands for the shaded region in Fig. |H] where the binomial 
probability is maximum.) The priors do not add up to unity yet; we shall comment on this 
point in a moment. 




FIG. 8: FRW cosmology: the values of (^mo> ^ao) m the shaded region maximize the binomial 
probability. 



c. Comparison 

We are now in a position to evaluate the final probability to affix to each model Mj 
as the "degree of belief" that, taking into account all the information in our possession, 
the hypothesis M i is true. All we have to do is to take each prior P(Mj) listed in Table 
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TABLE IV: Bayesian priors. 

Model Free parameters Prior 

Fixed point model (z tr > z max ) none 1/2 

Extended fixed point model 1 1/12 

FRW with n A0 = 1 1/12 

FRW with O A0 + n M0 = 1 1 1/12 

FRW with O,yio and SIao i n region 1 2 1/8 



IIV1 multiply it by the measured binomial likelihood P(DS\Mi) in Table ITTT1 calculate the 
normalization factor M = P{DS\M i )P(M i ) 1 and eventually compute the final probability 
P{Mi\DS) according to (jXTJ. 

Our results are summarized in Table Ivl The fixed point model is found to perform better 
than any standard cosmology. It has a posterior probability of 37.5%, while the best-fit FRW 
model achieves only 25%. Clearly these two numbers are not too different, but probably it 
is safe to say that the fixed point model is at least as likely to be the correct theory of the 
late Universe as is the best standard cosmology. Also it should be kept in mind that the 
Bayesian result always depends on the choice for the priors. While our choice seems natural 
and has been used already in similar applications [20] it is clear that others are conceivable 
as well. 

Please note that the numbers in the last column of Table IIVI do not add up to unity 



TABLE V: Final probabilities. 



Model 


Confidence model is true 
(after 92 SNe la) 


Fixed point model (ztr > z m ax) 


0.36 


Extended fixed point model 


0.16 


FRW with ft A0 = 


0.08 


FRW with fi A o + ^mo = 1 


0.16 


FRW with Omo an d Fl\o rn region 1 


0.24 



41 



yet. This is because we allow for models with more than 2 free parameters, whose prior 
probabilities, added up in the usual manner, ^2 N °^ 3 P(Mi), cooperatively make up for the 
missing 1/8. However, we can safely ignore these cosmologies in the calculation of the final 
probabilities, because they are strongly penalized by their small priors, l/2 Ni+l with iVj > 2, 
and their inclusion in the computation does not significantly alter the results. 

We can give a quantitative estimate of this assertion as follows: let us include in our 
discussion more and more cosmological models, Mj (i = 5, ■ ••), each characterized by iVj 
free parameters, with iVj > 2. Their priors P(Mj) are certainly such that 

P{Mi) < ^ (<> 4) (A4) 

On the other hand, their experimentally determined probabilities P(DS\Mj) are subject to 
an upper bound, too, because they cannot exceed the maximum of the binomial likelihood 
function P(n, 92), which is obtained for n = 46: 

P{DS\Mi) < P(46, 92) = 0.08 (A5) 

The increment in the normalization factor M = £\ P(DS\Mi)P(Mi), due to the new inclu- 
sions, is rather irrelevant: 

+oo +oo +oo . 

5N=Y, P{DS\M i )P{M i ) < 0.08 P{Mi) = 0.08 ^ ^vhT = 0.01 (A6) 

i=6 i=6 Ni=3 

The value of Af, for the first 5 models only, amounts to 0.042. Therefore, the effect of 
including all the possible cosmologies with N{ > 2 is a shift of the value of M from 0.042 to 
(at most) 0.052. 

By eq. (|A1[) . we also need to recompute the final probabilities P(Mi\DS). For the first 
5 models, we have: 

r-m»S) - P{DS J^l P J Mi) - pM< - M ^WTW < A7) 

Since 

< 5N < 0.01 (A8) 

the change in the final probabilities in Table IVI reduces to a simple rescaling by a factor that 
satisfies 

P^fMADS) , k , 

°' 81 ~ pAmads) £ 1 < A9 » 
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In particular, we observe that, irrespective of what the new models are (as long as N{ > 2, 
of course) or how well they compare to the data, the IRFP model will always have 

P ncw (Mi\DS) > 29% (A10) 

We might wonder how this figure compares to the new final probabilities for the other 
cosmologies. As for the four models with W < 2, since P old (M 1 |DS) > P old (Mi\DS), 
i = 2, 3, 4, 5, eq. (|A7J) shows that 

P ncw (M 1 | J D5) > P ncw (Mi\DS) (All) 

in all cases. (The final probabilities are merely scaled, so their ratios remain unchanged.) 

As a last step, we must make sure that some of the new models do not acquire a final 
probability which is larger than P ncw (Mi\DS). This possibility is easily ruled out by using 
(|A4|) . (HD and flUJ): 

r-WDS) = «^ < ^ = 0.096 (i = 6,..) (A12) 

Therefore, the new models cannot achieve a confidence level of more than ~ 10%. This 
upper bound is based entirely on the fact that they contain a number of parameters which is 
larger than 2. Our final figures, reported in Table El reflect precisely these considerations, 
which fully entitle us to ignore any model with Ni > 2. 
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